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Abstract 

We study a massless Dirac particle with PT symmetric non-Hermitian Rashba inter¬ 
action in the background of Dirac oscillator potential to show the PT phase transition 
in a (2-1-1) dimensional relativistic system analytically. PT phase transition occurs when 
strength of the (i) imaginary Rashba interaction or (ii) transverse magnetic held exceed 
their respective critical values. Small mass gap in the spectrum, consistent with other ap¬ 
proaches is generated as long as the system is in the unbroken phase. Relativistic Landau 
levels are constructed explicitly for such a system. 
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1 Introduction 


Consistent qnantum theory with real energy eigenvalnes, nnitary time evolution and 
probabilistic interpretation for combined Parity (P) and Time reversal (T) symmetric 
non-Hermitian theories in different Hilbert space equipped with positive dehnite inner 
product has been the subject of intrinsic research in frontier physics over the last one 
and half decades n-0 The huge success of the complex quantum theory has lead to 
extension to many other branches of physics and has found many applications |1]-|28]. 
Such non-Hermitian PT symmetric systems generally exhibit a phase transition 0-0 ( 
or more specihcally a PT breaking transition ) that separates two parametric regions (i) 
region of the unbroken PT symmetry in which the entire spectrum is real and eigenstates 
of the systems respect PT symmetry and (ii) a region of the broken PT symmetry in 
which the whole spectrum (or a part of it) appears as complex conjugate pairs and the 
eigenstates of the systems do not respect PT symmetry. The study of PT phase transi¬ 
tion has become extremely important due to the fact that such phase transition and its 
rich consequences are really observed in variety of physical systems [in]-[T6]. However 
most of the analytical studies on PT phase transition are restricted to one dimensional 
non relativistic systems. Few groups nn-iiHi have studied PT phase transition in higher 
dimensions for restricted class of non-Hermitian potentials, in non-relativistic situations. 
In the present work we consider a relativistic PT symmetric system in (2-1-1) dimensions 
to study the PT phase transition and its consequences. 

As a model we consider a relativistic system of a massless Dirac particle in (2-1-1) 
dimension [22]. Such systems have become extremely important due to the discovery of 
graphene, a single layer of the carbon atoms arranged in a honeycomb lattice. The study 
of graphene has created an avalanche in the frontier research over the past decade due to 
its exceptional electronic properties and potentially signihcant applications [30]-[35]. The 
low energy excitations of the atoms in the graphene obey massless Dirac equation in (2-1-1) 
dimensions. The two dimensional graphene sheet is essentially a zero gap semiconductor 
with Fermi level Ef precisely at i? = 0 with a linear dispersion relation given hj E = hvfk 
where Vf is the Fermi velocity and k is the carrier wave vector. However a small gap has 
been found experimentally in the graphene spectrum [36]. Such a mass gap is generated 
by breaking the inversion symmetry either (i) by considering spin-orbit interaction (e.g. 
Rashba interaction) or (ii) by the coupling the system with quantized circularly polarized 
helds. In this present work we consider a massless Dirac particle with PT symmetric 
imaginary Rashba interaction [37] in the background of Dirac oscillator (DO) potential 
[38]-[IS] to demonstrate the PT phase transition in (2-1-1) dimensional relativistic system. 
We analytically show that PT phase transition occurs in such a system when the strength 
of (i) the imaginary Rashba interaction or (ii) the transverse magnetic held exceeds a 
critical value. However it is challenging to make a realistic experimental setup to observe 
the PT-phase transition for the system of graphene with imaginary Rashba interaction. 
Small mass gap which is consistent with other approaches are also generated as long as 
the system is in PT unbroken phase. The mass gap vanishes at the phase transition 
point and becomes imaginary in the broken phase of the system. We construct the 
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solutions for the other valley of the system and show the PT phase transition in the 
similar fashion. Relativistic lowest order Landau levels are constructed explicitly for the 
system. Connection of this system to well-known Jaynes-Gumming (JC) |16]-[17] model 
in optics is also indicated. 

Now we present the plan of the paper. In section (2) we introduce the model and its 
symmetries. PT phase transition and generation of mass gap in the spectrum are discussed 
in section (3). Section (4) is kept for conclusions. 


2 The Model and PT Symmetry 

The massless Dirac particle with imaginary Rashba interaction in the transverse static^ f 
magnetic held in the background of DO potential is described by the Hamiltonian 

H = Vf[a.{fi — iKiffi) + x n).x, (1) 

where H = p and RTi, A are real constants and u/ is Fermi velocity. 

We consider the magnetic held along z-direction, B = B^k and choose the vector potential 
as H = (—0) in the symmetric gauge. The term linear in r is DO potential |38] as 
in the non-relativistic limit, it reduces to simple harmonic oscillator (SHO) potential with 
strong spin-orbit coupling. Initially, DO was introduced in the context of many body 
theory, mainly in connection with quark conhnement models in QCD [39]. Later the 
subsequent studies have revealed several exciting properties connected to the symmetries 
of the theory and DO has found many physical applications in various branches of physics 
M-HS]. Non-Hermitian version of Dirac equation is also considred earlier.In particular, 
(2 -|- 1) dimensional massless Dirac equation in the presence of complex vector potential 
and its bound state solutions are discussed in [5l]. Exact solutions for the bound states 
of a graphene Dirac electron in various magnetic held with translational symmetry are 
obtained in [55] . 

The Hamiltonian in component form is written as 

H = Vf{ajl^ +ayliy) - iKiVfiyOxX + ayy)(3 + - ayli^) (2) 

where, ^ and H^ = -t- 

Now we discuss the symmetry properties of this model. Parity transformation. 



is an improper Lorentz transformation (i.e. detH = —1) and is dehned in two alternative 
ways in two dimension as, 

Pi\x —> -X, y — >y, Pa; — -px, Py —t py. (3) 

^ Various solutions of Dirac equation in presence of inhomogeneous[45]-150) as well as time dependent 
[5T]-[S3] magnetic field have also been obtained for hermitian model 
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X, y 


-y, Px 


(4) 


P 2 -. X —> X, y —> -y, ^ p^, py —^ -py. 

Under such parity transformations Dirac wavefunctions transform as 


Pitlj{x,y,t) = aytlj{-x,y,t) 

P 2 'ip{,x, y, t) = a^ilj{x, -y, t) (5) 

such that free Dirac equation remains invariant under Pi and P 2 . The time reversal 
transformation (t — > —t, i — > —i, p^ — > —Px, Py — > ~Py) (2+1) dimension 
Dirac theory is dehned as T = iayK, where K is complex conjugation operation such 
that, TiIj^x, y, t) = iayK'il){x, y, —t) = iay'ijj*{x, y, —t). It is straight forward to check that 

+ aylly) - iKiVf{a^x + ayy)/3 - iX{a^Ily - ayll^) 

^ H 

This system, however is invariant under both PiT and P 2 T as 

PiTH{PiT)-^ = H 
P2TH{P2T)-^ = H 

A massless relativistic system has two degenerate but inequivalent solutions known as 
valleys of the system. For the present system in this work there exists another Dirac 
equation corresponding to the other valley. The Hamiltonian corresponding to other 
valley is obtained by taking a time reversal transformation of H [56l [57] as, 

H = THT~^ = Vf{axflx + (Xy^y) - iKiVf{axX + ayy)(3 - i\{axtiy - aytix) ( 6 ) 

where tlx = Px + and tly = Py — . It is straight forward to show that H W and 

[if, PiT] = [ii, P 2 T] = 0. Now in the next section we hnd the solution of Dirac equation 
for both the valleys to show the PT phase transition and mass gap generation in such 
systems. 


3 PT phase transition and generation of mass gap 


In this section we solve the Dirac equation with Hamiltonian H and H corresponding to 
the valleys for the massless particle explicitly to demonstrate PT phase transition in (2+1) 
dimensional relativistic system. In both cases we found that mass gap is generated due 
to imaginary Rashba interaction in the background of DO potential as long as the system 
is in the unbroken phase. We further show how the different solutions of Dirac equation 
for both the valleys in this particular model are interrelated. We begin with writing the 
Hamiltonians in Eqs. ([2]) and ([6]) in a compact form on a complex plane z = {x + iy) as. 


/ 0 AUz + iCiZ \ ~ 

V BU, + iC2Z 0 y) ’ ^ 


0 BUz + iC2Z 
AUz + iCiZ 0 


(7) 
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where, A = 2(vf — A), B = 2(vf + A), Ci = KiVf — {vf — A)^ and 
C 2 = —KiVf + {vf + A)^ are constants and canonical conjngate momenta in complex 
plane and are defined as = -ih-^ = \{px — Wy), = ^{px + ipy) 

with the properties 

[z, nj = ih] [z, n^] = 0; [z, nj = O; [rig, n^] = O; [z, n^] = ih (8) 


3.1 Solution for the system H 

Now we first present the solntions corresponding to the Hamiltonian H. To find the soln- 
tion of the Dirac eqnation corresponding to H we assnme a solntion of the corresponding 
Dirac eqnation B'l/^ = Etf) in the two components form as, 

(9) 

Dirac eqnation is then written in components form as 

{AYiz + iCiz){BYi2 + iC2z)(f) = E^(f) ( 10 ) 

(HH, + iC2z){An, + iC^z)x = E\ (11) 

Now we look for a solution of the kind 

X = v{z,z)e^^^^ (12) 

Substituting these in Eq. lfTOj) or in Eq. lfTTD and comparing the coefficient of {zz) in both 
sides we obtain di = ^ or For hi = ^ the Eqs. ffTOD and (fTTD reduce to, 

- ABh^^ + Kzp- -E^p = 0 (14) 

ozoz oz 

Similarly for di = the Eqs. ffTOj) and ffTTj) reduce to, 

- ABh^-^ - Kz§: - = 0 (15) 

ozoz oz 

- ABh^^ - Kz^ - {E^ + K)p = 0 (16) 

Ozoz Oz 

We find the solution of the Eqs. ffTSD and ffTTlh i.e. for di = -^ ) as 

= al^z^ 

Pn = with, El = {n + 1){AC2 - BCi) 
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where are real constants. Thus the general solution for the level for the case ^ 

is 




irjn 


= ai 


iz 


z"- 

n+1 


^ TP + 


(17) 


In exactly similar fashion we obtain the solution corresponding to di = as, 


i’i\z,z,t) = ai^ { lessee s 


^zz-iE„t 


, with, El = {n + l){BC, - AC 2 ) = -El (18) 


3.2 Solution corresponding to other valley 

To hnd the solution corresponding to the Dirac equation for H we observe that, 

AoB 

H - > H 

C 10 C 2 

Interchanging ^4 -H- 5, Ci -H- (72 is equivalent to A —)■ — A, i7i —)■ —Ki and Bq —)■ —Bq. 
Thus it is straight forward to obtain the solutions for the Dirac equation corresponding 
to H which describes their valley, by using these changes of parameters in the solutions 
corresponding to Dirac equation with Hamiltonian H. The solutions for the systems H 
are given as. 


~ C 

For = ^ , '4)l{z,zA) = a: 


/ I Z^~^^ \ ALzz —-E t 

iz"^ 


For di = — , ^l\z,z,t) = dl^ f .^n+i 


^2 _i_ fr i 


(19) 


( 20 ) 


where El = {n + 1){AC2 — BCi) = —El. We will be using these results to discuss PT 
phase transition and generation of mass gap in the next subsection. 


3.3 PT phase transition 


The energy eigenvalues for the solutions corresponding to di = di = ^ in both the valley 
are. 


En = ±\j {n + 1) 
This indicates a mass gap 


2{v) - A2) 


B^eh 


- AKivlh 


( 21 ) 


Aj = Aiv'j - X‘)!^ - TKWjn (22) 

is generated between the positive and negative energy solutions due to the interactions 
present in this theory. Now we proceed to show that this system passes through a PT 
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phase transition depending on the different parametric values. The energy eigenvalues 
are real for these solutions when 


(i) < tt(l - OR 

As long as En is real and are an eigenstate of both PiT and P 2 T 

PiTi’i(z,z,t) = ( 4i;i( A,-1) ) 

V / 


(23) 




Similarly, 


(24) 


PzTiii(z,z,t) = 


G ^%G y 


Z) \ ^ ^ 

-ir]n{z,z) I 


= —la^ 


-irjr 




^ ^^ZZ-iEr^t _ jJl 


irjn 




(25) 


Similarly we obtained for the other valley as, 

PiT'ipi{z,z,t) = -i{-l)^'tpi{z,z,t) 

P2T'ipi{z,z,t) = -^l{z,z,t) (26) 

Thus for di = di = -^ system corresponding to both the valleys are in unbroken PT 
phase as long as the coupling for Rashba interaction is less or equal to Ac or strength 
strength of the magnetic held is greater or equal to a critical value Bq. Now we observe 
in this situation the other two solutions {ipn,'ipn) (i.e for di = di = ^) correspond to 
the broken PT phase as the eigenvalues for both the valleys are imaginary, En = EiEn 
with real En- In this case ijjn and tljn are not eigenstates of either PiT or P 2 T i.e., 

PGi^i^ ^ for ^ = 1,2 (27) 

as energy eigenvalue En is imaginary. Hence PT symmetry is broken spontaneously for 
di = di = ■^. Thus the system with di = di = ^ passes through PiT and P 2 T phase 
transition where either (a) strength of the imaginary Rashba interaction is equal to a 
critical value (Ac) or (b) the external magnetic held Bq is equal to a critical value {Bq). 

On the other hand if A > Ac, for hxed Bq or Bq < Bq for hxed A, the solutions for both 
the valleys for di = di = i.e. 'ijjpp,ippp are in the unbroken phase and the solutions 
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corresponding to di = di = i.e. are in broken phase. PT phase transition for 

di = hi = ^ is demonstrated graphically in Fig. 1 by plotting vs A for fixed value of 
Bq in Fig. 1(a) and by plotting En vs Bq for fixed A in Fig 1(b). 


E„ 



E„ 



(a) (b) 

Fig. 1: PT phase transition for the case for di = di = ^. Fig (a) Real (Solid 
lines) and imaginary (broken lines) parts of En have been plotted with Rashba coupling 
strength for fixed values of Bq = 100, Ki = 0.02, Vf = 0.01c, e = 1, h = 1, c = 137 for 
n = 0, 1,2, 3,4. Phase transition occurs at Ac = 1.33193. Fig (b) Similarly En is plotted 
with strength of magnetic field {Bq) for fixed values of X = 0.5, keeping other parameters 
same as in Fig 1 (a) for n = 0,1, 2, 3,4. Phase transition occurs at Bq = 6.32209. 

Now we would like to concentrate on mass gap Aq when the system passes a PT phase 
transition; 

(i) For di = di = ^, when A < Ac or i?o > Bq we see from Eq. fl22l) mass gap is 
always positive and becomes imaginary when A > Ac or Sq < Bq. Thus a mass gap is 
generated as long as system is unbroken PT phase and given by the expression in in Eq. 
fl22|) . Further in unbroken phase Aq varies as \/Bq for large enough magnetic held. This 
result is consistent with other approaches. 

(ii) For di = di = ^ when A > Ac or i?o < Bq the system is in unbroken phase and 
again real mass gap is generated which is consistent with other approaches. Thus mass 
gap is always generated in this theory which is consist with other approaches. In Fig. 2 
variation of mass gap is shown with respected to the Rashba coupling strength as well as 
strength of the magnetic held. Real mass gap for the massless systems is generated as 
long as the system is in PT unbroken phase. 







Ao Aq 




(a) (b) 

Fig. 2: Variation of mass gap with A and Bq for the case di = di = ^. Fig. (a) 
Real (Solid line) and imaginary (broken line) part of mass gap is plotted with strength of 
the Rashba interaction for fixed values of Bq = 100, Ki = 0.02, Vf = 0.01c, e=l, h = 1, 
A < Ac = 1.33193, mass gap is real and system is unbroken FT phase. Fig. (b) Real 
(Solid line) and imaginary (broken line) part of mass gap is plotted with the strength of 
magnetic field by keeping the strength of Rashba interaction fixed for a value of X = 0.5, 
Ki = 0.02, Vf = 0.01c, e=l, h = 1. Bq > Bq = 6.32209, system is in unbroken phase, 
and real mass gap varies as Aq oc y/^. 


Special cases : 

(i) In absence of DO {Ki = 0), system passes to broken phase when strength of 
Rashba interaction exceeds Fermi velocity X'^ > vj. The mass gap in that case (unbroken 

PT phase) is given by Aq = ^2(r;| — which is real as A^ < vj for the unbroken 

phase. Mass gap vanishes at the transition point A^ = vj. 

(ii) In absence of Rashba interaction (A = 0), the system is Hermitian, but can have 
real energy eigenvalues only when the magnetic held is sufficiently large i.e. Bq > 

the mass gap generated in this case is Aq = yj2v‘fh{^^ — 21^ which is always real for 

-®o > both the cases for large enough magnetic held Aq oc y/B^, which is consis¬ 

tent with other approaches. Again in this case, mass gap vanishes at the transition point, 
Bq = This critical value of magnetic held is same as obtained in jUj. The factor of 
2 is due to the choice of vector potential in the symmetric gauge. This discussion in valid 
for both the valleys. 
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Lowest Landau Levels (LLL) : 


The LLLs are obtained by putting the condition QiXo = 0, which implies 


(v4n^ + iCiz)xo — 0. 


We substitute Xo = z)e 


= 

. ^ zz. 
Ah 


to solve Eq. fl28ll and obtain 


( 28 ) 


dua{z,z) ^ 
dz 


(29) 


as the dehning rule for LLL in this system. We obtain the LLL in the coordinate repre¬ 
sentation, which is inhnitely degenerate as 

Xq{z,z) = ^ / = 0,1,2,.cx) (30) 


The monomials z^ with I = 0,1,2, ...oo serve as a linearly independent basis. The first 
and higher excited states in the coordinate representation are then obtained by applying 
Q\ given in Eq. (l30|) to the LLL repeatedly. The LLL corresponding to the other valley 
is obtained by putting (SII^ -|- iC 2 z)xQ = 0. Following similar calculation we obtain LLL 
Xo(z, z) = / = 0,1, 2, ...CX) and other higher excited states for this valley. 


Connection with JC Model ; 


We would like to point out the connection of this model with generalized version of 
well known JC model for massless particles in optics. The Hamiltonian corresponding to 
the simple version of (JC) model is given as 

Hjc = g{a+a + a_o)) + (31) 

where a± = iay) are usual spin raising and lowering operators. This model is 

widely used in optics to study the atomic transition in two level systems. 

Now the Hamiltonian H in Eq.(I7])can be expressed as. 


H = ^ 


0 Qi \ 

Qi 0 j 


(32) 


where, Qi = ^(HH^ J- iCiz) and QI = + iC 2 z) with K = {AC 2 — BCi)h. It 


is straight forward to check. 


Qi, QI 


= 1. Note Qi = Q 2 for real Rashba interaction. 
The Hamiltonian is not Hermitian but invariant under PiT and P 2 T both and is further 
written as, 

H = VK{a+Qi + (33) 


This is the extended version of JC model for particle with rest mass zero. It is very 
exciting to note how two completely different theories are related in this fashion. This 
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provides an alternative approach to study any massless Dirac particle (such as graphene) 
with Rashba interaction in the background of DO potential in the transverse magnetic 
held using JC model in quantum optics. For a real Rashba interaction, Qi = Q 2 and 
the extended JC model in Eq. fl33l) exactly coincide with JC model for a massless particle. 
Similar discussion is also valid for the other valley. 


4 Conclusions 

We have studied PT phase transition in a (2+1) dimensional relativistic system to enlight 
different aspects of this important phenomenon. For this purpose we have considered the 
massless Dirac particle with imaginary Rashba interaction in the background of DO poten¬ 
tial in (2+1) dimension. Such a system is non-Hermitian due to the imaginary coupling. 
We have constructed two different Parity- Time reversal transformations {PiT and P 2 T 
) in (2+1) dimensional relativistic system such that the system is invariant under both 
PiT and P 2 T. We have constructed different possible solutions corresponding to both the 
valleys for this system. The energy eigenvalues are real as long as the eigenstates respect 
PiT and P 2 T symmetry. We further investigate PiT and P 2 T phase transitions in this 
system to show that system passes from unbroken PiT and P 2 T phases to broken phases 
depending on the strength of the Rashba interaction or transverse magnetic held. The 
inter-relation of different solutions for both the valleys with respect to PT-phase transi¬ 
tion is presented. Even though graphene is zero gap semiconductor with linear dispersion 
relation, a small mass gap in its spectrum is observed experimentally. In our formulation 
small mass gap consistent with other approaches is generated for the system of massless 
Dirac particle when the system is in unbroken phase. The mass gap vanishes at the point 
of PT phase transition. We would like to point out that results on PT-phase transition 
presented in this paper is valid for any system with massless Dirac particle in presence 
of imaginary Rashba interaction in the background of Dirac oscillator. In particular our 
model has some relevance to the possible physics in graphene when the solution for both 
the valleys are combined. The lowest order Landau levels for this system have been cal¬ 
culated explicitly to show the inhnite degeneracy in such systems. The connection of 
such system with JC model has also been discussed. It will be exciting to correlate other 
exciting properties of graphene with PT phase transition. 
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